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ABSTRACT 


The  Monte  Carlo  fast-neutron  air-scattering  data 
presented  in  PZK-9-147,  Volumes  I  and  II,  have  been  inte¬ 
grated  to  obtain  the  angular  distributions  and  energy- 
spectra  for  a  point  isotropic  source  emitting  one  neutron 
per  second  at  a  given  energy  E^.  These  calculations  were 
performed  for  source -detector  separations  of  10,  35,  6ij., 
and  100  feet  and  for  initial  neutron  energies  of  0.33, 
1.1,  2.7,  U.O,  6.0,  8.0,  10.9,  and  1^.0  Mev. 
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1*  INTRODTJCTTOH 


All*  scattering  Is  the  principal  process  by  vhlch 
neutrons  leaving  e  sovirco  in  air  are  transported  to  other 
positions  some  diataaca  away*  In  order  to  make  comprehen¬ 
sive  shield-design  studies,  ona  must  know  the  energy  and 
angular  distributions  of  tha  scattered  neutron  flux  at  the 
positions  of  interest.  It  is  possible  to  design  and  apply 
Monte  Carlo  procedures  which  solve  this  problem  directly 
for  specific  cases;  however,  it  Is  often  more  practical  to 
use  the  Monte  Carlo  procedures  to  generate  data  in  a  suita¬ 
ble  parametric  forift  so  as  to  be  able  to  apply  the  results 
without  having  to  run  a  new  problem  for  each  new  source 
term. 

This  report  presents  parametric  air-scattering  data 
for  isotropic  sources  of  nonoenergetlc  fast  neutrons.  For 
source  energies  of  0*33#  l.l»  2,7,  i^.O,  6.0,  8.0,  10.9, 
and  lU«0  Mov,  the  angular  distributions  of  the  scattered 
flux  and  dose  rate  and  the  energy  spectra  of  the  scattered 
flvix  are  given  for  source-receiver  separation  distances  of 
13#  35#  64-#  *^hd  100  feet* 

The  data  presented  here  were  computed  from  the  results 
of  Wells*  Monte  Carlo  parameter  study  of  neutron  scattering 
for  directional  point  sources  in  an  infinite,  homogeneous 
fiiedlum  of  air  (Befs*  1  and  2)*  Ji  straight-forward  integration 
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procedure  was  formulated  and  programmed  in  Fortran  for  the 
IBM-70i|.  in,  order  to  carry  out  the  calculations.  The  pro* 
cedvire  was  designed  to  utilize  directly  the  pvinched  card 
output  of  Wells*  Monte  Carlo  calculations# 

Section  II  of  this  report  deals  with  the  geometry  of 
the  neutron-scattering  problem  while  Section  III  describes 
the  Monte  Carlo  data  \ised  In  these  calculations.  The  inte¬ 
gration  scheme  and  Fortran  procedure  are  discussed  in  Section 
IV  and  the  results  are  presented  in  Section  V  in  the  form 
of  tables  and  graphs  giving  the  angular  distributions  of  the 
scattered  flux  and  dose  rate  and  the  energy  spectra  of  the 
scattered  flux# 


II.  NEUTRON  SOATTERINa  GEOMETRY 


The  geometry  (Pig.  1)  consists  of  a  neutron  source  S 
located  In  an  Infinite,  homogeneous  medium  of  air*  A 
detector  D  of  unit  cross  section  is  located  a  distance  a 
from  the  source  3«  The  neutron  current  leaving  the  source 
is  described  by  a  polar  angle  K.  and  an  azimuthal  angle 
The  polar  angle  is  measured  with  respect  to  the  so\irce- 
detector  axis  while  the  azimuthal  angle  is  the  angle  between 
the  positive  y  axis  and  the  projection  of  the  neutron 
direction  on  the  y,z  piano*  Tho  azimuthal  angle  is  measured 
In  a  clockwise  direction* 

The  detector  angle  p  is  the  polar  angle  between  the 
direction  of  the  Incoming  neutron  flux  and  the  source^detector 
axis*  The  azimuthal  angle  at  the  detector  Is  defined  In 
the  same  manner  as  0*  The  neutron  current  of  energy  Ep 
moving  In  the  direction  at  a  point  q  on  the  surface 

of  a  unit  sphere  about  S  is  defined  as  5(K,^,Ep)*  The  num¬ 
ber  of  neutrons  passing  through  a  surface  element  dA  at  ^ 
per  unit  time  is  given  by  S  (K^j^jE^jdA. 
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flOURB  1*  SCAtTEHING  GEOMETRY 
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III.  MONTS  CARLO  7)ASk 

Since  the  generation  cf  the  Monte  Carlo  data  used  In 
tU©a«  calculations  has  been  thoroughly  described  by  VTells 
in  BeferenoQ  1,  only  a  brief  review  is  needed  here. 

The  cf facts  of  elastic  scattering,  Inelastic  scattering 
end  absorption  vere  taken  into  account  in  the  Monte  Carlo 
air -scattering  calculations.  The  neutron  cross  sections  for 
nitrogen  and  oxygen  used  in  these  calculations  were  taken 
fron  data  giving  by  Lustlg,  Goldstein,  and  Kalos  (Refs.  3 
and  I|.)  •  The  neutron  cross  sections  for  air  were  computed 
on  the  basis  of  a  mixture  of  nitrogen  and  22^  oxygen 
and  a  density  of  5*37  x.  10^^  atoms/co  or  0.1293  x  10”^  gra/cc. 
The  flux-to-dose  conversion  factors  F(e)  used  In  cotaputlng 
the  tissue  dose  rates  were  those  calculated  by  Burst  and 
Ritchie  (Fef.  5)* 

The  ealoxxlated  results  presented  in  Reference  2 
represent  the  angular  distributions  of  the  neutron  flux, 
the  angular  distributions  of  the  tissue  dose  rate ,  and  the 
energy  distribution  of  the  neutron  flux  at  the  detector  D 
for  a  point  monodlreetlonal  eouroa  8  emitting  one  neutron 
per  second  of  energy  in  the  given  direction  (K,^).  The 
calculations  were  performed  for  four  separation  distances 
of  10,  35,  6^,  and  100  feet,  eight  Initial  neutron  energies 
of  0.33,  1.1,  2.7,  6.0,  8.0,  10.9  and  1^.0  Mev,  and 
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«l^t  source  angles,  K,  of  l$f  30,  90,  120,  150,  and 

l80  degrees.  ^Chls  represents  a«. total  of  6I|.  different  source 
descriptions  and  four  source -detector  separations  which  nay 
be  combined  to  describe  any  arbitrary  source  distribution* 
The  scattered  flux  reaching  the  detector  is  sorted 
into  eighteen  equal  angle  intervals,  and  ten  arbitrary 
energy  intervals,  E^*  N(K,jZf,Ejj,a,p^,Bj)  is  defined  as  the 
fraction  of  neutrons  with  Initial  energy  Bq  and  direction 
paaalng  through  the  Incremental  area  dA  per  unit  time 
which  reach  the  detector  at  a  distance  a  in  the  detector, 
angle  increment  with  energy  in  the  energy  group  Ej* 
D(K,)2^,EQ,a,Pj^)  is  defined  as  the  neutron  dose  rate  due  to 
neutrons  with  initial  energy  and  direction  (K,)^)  which 
arrive  at  the  detector  D  a  distance  a  from  the  source  In 
the  detector  angle  inorenent  Since  there  le  symmetry 

about  the  source -detector  axis  and  both  source  and  detector 
are  located  in  the  infinite  homogeneous  medium,  angular 
distributions  at  the  detector  are  Independent  of  both  source 
and  detector  aalmuthal  angle*  Therefore,  the  data  are  not 
sorted  with  respect  to  this  angle,  and  the  variable  0  is 
dropped  from  the  expressions* 

The  following  quantities  are  tabulated  In  Reference  2s 
N(K,EQ,a,Pj_,Ej),  lT(K,Ep,a,EjJ,  N(K,Eo,a,p^) ,  D(,Bo,a,p^), 

BCK,EQ,a),  and  I)CK,Ejj,a) .  lI(K,EQ,a,Pj^,E j)  is  tabulated  in 
the  form  of  an  18  x  10  matrix  for  each  of  the  four  values 


of  a,  eight  values  of  E^,  and  eight  values  of  K«  In  addition 
to  the  l8  X  10  matrices,  th.e  energy  distribution  of  the  air- 
scattered  flux  is  given  by 


N(K,E^,a,Ej) 


(la) 


the  angular  distributions  of  the  air-scattered  flux  is  given 

hy 


10 

N(K,E^,.a,0^,E.)  .  C2a) 

j=*l 

Tho  do3»  rate  Is  calculated  and  stored  as  a 

function  of  K,a^,a  and  0^.  It  la  printed  as  a  function  of 
angle  0  for  each  combination  of  and  a. 

The  total  scattered  flux  and  dose  rata  are  given, 
respectively,  by 


N<K,Eo,a,pi) 


r 


N(K,EQ,a) 


and 
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D{K,E^,»)  =  ^  P(K,Bq,a,0j^)  .  (5a) 

The  Units  Of  the  scattered  flxixes  are  neutrons/on^-sac  per 
source  neutron/sec  and  those  of  tho  scattered  dose  rates  are 


rom/lur  per  sotarce  neutron/sec. 


The  total  scattereci  flux,  XrOiti  %  source  SCKjE^)  teaching 
tjie  (detector  D  at  a  distance  a  Is 

X  211 

H(B^,a)  *=  /  /  SCK,Ej^jN(K,B^,a)  siniC  dK  (6a) 


IT 

2^  f  S(K,BQ)lT(K,Eo,a)  elnK  dE. 


Stmllatly,  the  total  scattered  dose  rate  Is 


DCB^,a)  ”  /  /  S(K,BjD(K,Ejj,a)  sinK  AK  Ai 


»  2%  j  S{K,S-)D(K,Eo,a)  sinK  dK. 


Slallarly*  the  angular  distributions  of  flux  and  dose  rate 
are  given,  respectively,  by 


lI(So,a, 


IT  2n 


"7/ 


5(K,E->lI(E,Eo,a,pJ  sinK  dK  d^  (8a) 


23C  /  S(K,Eo)II(K,E-,a,p)  sinK  dK,  (8b) 


%  2% 


■//  S(K,EQ)D(K,B^,a,(3)  sinK  dK  (^a) 

ft 

=  2ft  J  3 (E,Eo)  D  (r,EQ , ft ,  P)  sln!i  dE.  (Vo) 
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The  energy  distribution  Of  the  scattered  fltuc  is  given  by 


NCa,EQ‘E) 


SCK,B^3U(K,Eq 

CK,Bo3ir(K,E^, 


,a^S}  slnS  dK 
a,E3  sinK  dK« 


(10a) 


(lOb) 


An  analysis  of  the  Monte  Carlo  data  shovad  that  the 
standard  deviation  of  the  total  scattered  fltuc  iKKjE^ja)  was 
less  than  0.1  N(K,E^,a)  in  60t2^  of  the  problems,  less  than 
0.15  KCKjE^ja)  in  91*4-^  of  the  problems,  and  in  only  Z,73% 
of  the  problems  did  the  standard  deviation  exceed  0*2 


4 


IV.  INTSaHATlON  OP  THE  HONTE  CARLO  DATA 

SiriGQ  the  quantities  llCK,EQ,a.p) ,  D(K,E^^,a,,p) ,  and 
NCK,EQ,a,E)  are  not  given  In  terms  of  analytic  functions, 
it  is  necessary  to  devise  some  integration  scheme  that  Vill 
produce  accurate  results*  Eor  the  purpose  cf  numerioal 
integrations.  Equations  6h,  7h,  8b,  9t,  and  10b  become.  In 
turn 

II(E^,a)  a  S(K^,Ej^)N(K^,E^,a)  sinK^  (AK^),  (6c) 

6 

D(EQ,a)  a23:\  S )D  (K^,E^,a)  sinK^  (AK^),  (7c) 

n=»r 

_8, 

»(Eq,II,Pj^)  "  ^  S(K^,Eo)N{K^,E^,a,Pj)  sinK.^  (AK^),  (8c) 

PCB^,,a,Pj^)  a  211  2  S(K^,B^)DCK^,B^,a,Pj)  slnK„  (AK^),  (9c) 

n* 

N(E^,a,Ej)  =  211  S  S(Efj,E^)NCK^,B^,«,Bj)  fiial^  (AK^) .  (lOc) 
vhe  re  1^1 ,  2,  ••••  ,  l3  and  j  2,  ••*  ,  XO  * 
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Th©  calculations  presented  here  Were  performed  for  a 
point  isotropic  source  emitting  1  neutron/sec  with  energy  Eq» 
Therefore , 

S(K,Eq)  *  constant  =  (11) 

The  values  of  N  and  D  were  taken  from  the  tables  in  Refer¬ 
ence  2* 

The  suffiinatlons  In  Equations  6o  and  7o  were  performed 
using  a  ttiachlne  code  designed  primarily  for  matrix  calculations 
but  adaptable  for  those  oalculations* 

An  IBM  Fortran  procedure  was  written  to  provide  for  the 
integration  of  E<iuatlona  8o,  9o»  ahd  lOo  with  a  noalsotroplc 
soTxree  tern  depending  only  on  the  polar  angle  K*  The  Isotropic 
case  presented  here  is  a  special  case  of  this  more  general 
source  term,  A  copy  of  the  Fortran  program  is  shown  in  the 
Appendix.  .  Certain  subroutines,  such  as  SETUP,  END 
LIB  1,  LIB,  and  END(2,1),  are  special  (Jeneral  Dynamlcs/Fort 
Worth  subroutines.  However,  the  basio  program  should  remain 
the  same  anywhere. 

The  values  of  used  In  the  numerical  integration 

ware  chosen  after  considering  which  histogram  would  best 
represent  the  smooth  curve  for  Integration  purposes*  The 
histogram-fit  tends  to  ovsrestlmate  in  the  first  20  degrees 
•nd  underestimate  In  the  last  20  degrees* 


* 

The  Validity  of  the  Ititegration  scheme  was  checked 
by  comparing  the  results  of  htimerical  Integration  vlth  those 
obtained  using  a  planimeter  on  smooth  curves  drawn  through 
the  eight  points#  In  the  seven  cases  chosen  at  random  for 
comparison,  the  numerical  integration  overestimated  the 
result  by  less  than  In  all  oases,  the  numerical  inte¬ 

gration  overestimated  the  result. 


V*  hesults 


Eg,tiatsions  6o,  7o,  9o,  and  10c  liav©  baan  numeJ^ically 

Integrated  fon  a  point  isotropic  sotiroe,  The  results  are 
presented  In  this  section  in  both  tabulated  and  graphical 
form , 


The  scattered  neutron  flux,  lI(B^,a,Ej),  Tables  II 

through  V,  and  the  neutron  flux  per  Mev, 

Table’s  VI  through  IX,  have  been  tabulated  as  functions  of 

Bjj,  the  Initial  energy,  and  ..  The  subscript  n  indl- 

i  J 

cates  that  CE^)j  is  the  lower  limit  of  the  energy  interval* 
Ce^)j  for  each  Eq  la  shown  In  Table  I,  The  minlwurii  energy 
cutoff  for  each  source  energy  is  given  In  Table  I  by  (e^j)  j 
for  each  source  energy.  AE^  Is  the  width  of  the  energy 
group.  The  scattered  neutron  flux,  Is  reported 

in  units  of 

neutrons/ cm^»sec 
source"  neutron/ sec 

and  the  neutron  flux  par  Mev,  h  CEo,a,Bj)/AEj  in  units  of 

neutrons/cm^-sec/jtev 


source  newtron/sec 


The  angular  distributions  of  the  scattered  neutron  flux 
end  dose  rate  in  a  10®  interval  of  and 

are  tabulated  as  functions  of  end  in 
Tables  X  through  XIII  and  XIV  through  XVII,  respeotivoly, 


m 

*  ii 

Vher©  is  the  iippej?  limit  of  the  angulaj*  interval* 

Ihe  units  here  are 

neutron a/cm^-se a  in  10*^  interval 
source  xieutron/seo 


and 


rem/hr  In  10^  Interval 
source  neutron/sea 


The  quantities  NCE  ,a,S,)/AE.|  UCE  .a,p  ),  I>CE 
N(S^,a) ,  and  are  shown  in  Figures  2-9#  10-17#  18-25# 

26  and  27i  respectively. 

The  angular  distributions  are  plotted  as  funotion  of 
the  midpoint  cf  the  angular  interval  to 

In  Flguroa  2  through  9#  the  energy  spectrum  is  plotted 
against  the  midpoint  of  the  energy  group  E^, 

The  total  scattered  flux  and  dose  rate*  NCB  .a)  and 

0 

D(Eo,a),  respectively#  are  shown  in  Figures  26  and  27  as 
functions  of  E^^  for  each  of  the  four  separation  distances 
considered* 

The  scattered  dose  rates  for  source  energies  of  1,1  Mev 
or  greater  result  from  neutrons  with  energies  greater  than 
the  minimum  energy  used  for  each  source  energy,  but  the 
scattered  doae  rates  for  the  0*33  Mev  source  are  those  re¬ 
sulting  from  neutrons  with  energies  greater  than  0*22  Mev, 
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TABLE  XIII  .  ANGULAR  DISTRIBUTION  OF  TOTAL  SCATTERED  NEUTRON  FLUX 

Separation  Distance  -  100  Feet 


Source  Energy  (Mev) 
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TABLE  XV  ANGULAR  CISTRIBUTION  OP  TOTAL  SCATTERED  NEUTRON  DOSE  RATE 

Separation  Distance  -  35  Feet 
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TABLE  XVI  •  ANGULAR  DISTRIBUTION  OF  TOTAL  SCATTERED  NEUTRON  DOSE  RATE 

Separation  Distance  -  64  Feet 
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FIGURE  2.  ENERGY  SPECTRA  OF  TOTAL  SCATTERED  NEUTRON  FLUX 

Initial  Energy  0.33  Mev 
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Scattered  Neutron  Flux  [(neutrons/cm2-sec-Mev)/( 
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FIGURE  4.  ENERGY  SPECTRA  OF  TOTAL  SCATTERED  NEUTRON  FLUX 

Initial  Energy  2,7  Mev 
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FIGURE  5.  ENERGY  SPECTRA  OF  TOTAL  SCATTERED  NEUTRON  FLUX 

Initial  Energy  4.0  Mev 
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Scattered  Neutron  Flux  ^neutron3/cm^-3ec-Mev)/(30urce  neutron/3ec 
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FIGURE  6.  ENERGY  SPECTRA  OP  TOTAL  SCATTERED  NEUTRON  FLUX 

Initial  Energy  6.0  Mev 
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Scattered  Neutron  Flux  [(neutrons/cra^-sec-Mev)/(source  neutron/sec 
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FIGURE  7.  ENERGY  SPECTRA  OP  TOTAL  SCATTERED  NEUTRON  FLUX 

Initial  Energy  8.0  Mev 
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Scattered  Neutron  Flux  [( 


FIGURE  8.  ENERGY  SPECTRA  OF  TOTAL  SCATTERED  NEUTRON  FLUX 

Initial  Energy  10.9  Mev 
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FIGURE  9.  ENERGY  SPECTRA  OF  TOTAL  SCATTERED  NEUTRON  FLUX 

Initial  Energy  l4.0  Mev 
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FIGURE  11.  TOTAL  SCATTERED  NEUTRON  FLUX  VS.  DETECTOR  ANGLE 

Initial  Energy  1.1  Mev 
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Scattered  Neutron  Flux  [(neutrons/cm^-sec-lO'^  Interval )/( source 
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FIGURE  12.  TOTAL  SCATTERED  NEUTRON  FLUX  VS.  DETECTOR  ANGLE 

Initial  Energy  2.7  Mev 
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Scattered  Neutron  Flux  f(neutron3/cm2-3ec-10O  Interval) /Csource  neiitror^aaejj 
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?IGUEE;  14.  TOTAL  SCATTERED  NEUTRON  PLUS  VS.  DETECTOR  ANGLE 

Initial  Energy  6,0  Mev 
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FIGURE  15.  TOTAL  SCATTERED  NEUTRON  FLUX  VS.  DETECTOR  ANGLE 

Initial  Energy  8.0  Mev 
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Scattered  Neutron  Flux f(neutrons/cm2-3ec-10O  interval)/( source  neutron/sec^ 
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FIGURE  16.  TOTAL  SCATTERED  NEUTRON  FLUX  VS.  DETECTOR  ANGLE 

Initial  Energy  10.9  Mev 
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FIGURE  17.  TOTAL  SCATTERED  NEUTRON  FLUX  VS.  DETECTOR  A; 

Initial  Energy  14.0  Mev 


5S 


Scattered  Neutron  Dose  Rate  Prem/hr-10®  interval) /(source  neutron/sec 


Detector  Angle  /3  (degrees) 


FIGUHS  l8.  TOTAL  SCATTERED  NEUTRON  DOSE  RATE  VS.  DETECTOR  ANGLE 

Initial  Energy '0.33  Mev 


60 


Detector  Angle  (degrees) 


FIGURE  19.  TOTAL  SCATTERED  NEUTRON  DOSE  RATE  VS.  DETECTOR  ANGLE 

Initial  Energy  1.1  Mev 
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FIGURE  20. 


50  75  lOO  125  150  175 

Detector  /Vngle  yS  (degrees) 

TOTAL  SCATTERED  NEUTRON  DOSE  RAT®  VS.  0  .'ECT''."  ../.ILK 
Initial  Energy  2.7  Me.'-s 


62 


Scattered  Neutron  Dose  Rate  [treiVhr-10‘^  lnt3i’Val)^( 
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FIGURE  22,  TOTAL  SCATTERED  NEUTRON  DOSE  RATE  VS.  DETECTOR  ANGLE 

Initial  Energy  6.0  Mev 
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Scattered  Neatron  Dose  Rate  ({ rem/hr- 10°  interval VCtOUrce  netftron/sec^ 


Scattered  Neutron  Dose  Rate  j(rem,/hr-10°  interval)/( source  neutron/set: 


Til 


50  75  100  125 

Detector  Angle  (degrees) 


PTOURE  25.  TOTAL  SCATTERED  NEUTRON  DOSE  RATE  V3.  DETECTOR  ANGLE 

Initial  Energy  l4.0  Mev 
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Scattered  Neutron  Flux  [{ neutrons/ 


FIGURE  26,  TOTAL  SCATTERED  NEUTRON  FLUX  VS.  INITIAL  ENERGY 


68 


source  neutron/sec 


APPENDIX 

J?-55  POBTRAN  STATEMENTS 


71 


APPENDIX 


'55  FORTRAN  STATEMENTS 

7171  BEG  XTR«  R5S  JACK.  GRIGSBY  CCfiRO.L  DIFFEV— 3987J 
IZE  EZEC1  EZEI:i2  EZeD3  EOF  EZfOS  EZ£B1flE2EB3  EZEB4 
;  IBM  704  PROCEDURE  R55 

;  IHTEGRflTIOH  OF  MONtE  CrRLO  „,-;i.rULflT IOMs  OF  FfIST 

'  IW  fllR  FOR  HOM- ISOTROPIC  SOURCES 

DIMENSION  ZCB:>.  DZCEO,  FLUKfi  (:576,85.  DOSEfI 

1  FLUXe  C3-20,8j,  SC8>8>>  SRZCSI'- 

2  SS(:576,85  ,  FDC57B5.  SOS  <:iS-4,3;:«» 


FflO  C57S,! 


flibcs'j.  RZCSJ 


COMMON  Z  ,  DZ  ,  FLUXO  ,  DOSEfi  ^  FLUXE  »  S 
1  fl  ,  SS  .•  FD  ,  SOS  .  FPEi  f  FLIB  ,  RZ 

IQ  CALL  Liei  CM? 


GO  TO  Cl 

1  ,50?  , 

M 

It 

call  lib 

C3HR55, 

ED> 

read  500 

,  FLIP 

IF  C  FL 

IP  ?  14. 

*  ** 

1  w 

13 

BACKSF'AC 

E  9 

G  0  T  0  5  0 

1  4 

FLIP  = 

MODF  C 

ED  . 

»0.  0 

LIBT  = 

FL  I P 

GO  TO  Cl 

5,20.25.. 

30  f  3 

5?  , 

LIBT 

40 

FORMAT  C 

6E 1  0 . 4? 

1  5 

READ  40. 

C  ZCJ? 

,  j= 

J.S  > 

XID1  =  ED 

GO  TO  1 1 

20 

R  E  R  D  4  iD  t 

CDZCJ 

-1 

=  1 , 8 

XID2  = 

EO 

GO  TO  11 

25 

READ  40, 

CCFLLIXA  Cl 

>  J  j*  » 

j=i . 

8>  1  = 

XID3  = 

ED 

GO  TO  11 

30 

READ  40, 

CCDOSEA  Cl 

J  * 

j=i . 

II 

CO 

XI04  = 

ED 

GO  TO  11 

35 

READ  40. 

CCFLLIX 

IE  Cl 

.  ..1?  , 

j=i . 

8?  I  = 

XID5  = 

ED 

GO  TO  11 

50 

CALL  SETUP  C3HR5 

5, XID? 

READ  500 

.  CFLIE! 

CN?  . 

N=1  . 

5i> 

iOO 

FORMAT  C 

5F7.0> 

IF  OilDI 

-  FLIEC 

1?? 

55 » 

51  . 

IT  IT 
sJ 

51 

IF  CXID2 

-  FLIBC 

2';:) 

55  * 

C  C 

52 

IF  CXID3 

-  FLleC3?? 

55.' 

53. 

e-  C 

V.'  -k' 

53 

IF  CXID4 

-  FLIEIC4?? 

55  r 

54, 

54  IF  CXID5  -  FLIBC5??  55.  58.  55 

55  PRINT  56 

56  FORMAT  C54H  LIBRPRV  DECKS  CALLED  FOR  IN  PROBLEM 

CALL  EN09  .  . . . 

58  READ  59,  CN0P1  ,  N0P2  ,  N0P3  > 

59  FORMAT  C3I3? 

IF  CN0P1  +  N0P2.  +  N0P3  ?  1  66o.  1  00pj._  6 0 _ 

60  READ  40.  CCSCI  .J),  J  =  1 , 8? ,  I  =  1 , 8  ? 

DO  70  J=1 . 8 

RZCJ?  =  0.01745329  ♦  .  ZCJ? 

SR2CJ?=  SINF  CRZCJ?? 

ACJ?  =  SRZCJ?  *  DZCJ? 

70  CONTINUE  _  .  _ _ 


0001  ,R5,5 

EZEB1 AEZEB5 

0002  R55' 

neutron  SCAT.T.ER.I-NG0003  .'R,5.5 
i?004  R55 
0005  R55 

C576,,  S?.,  .  0006.  R55 

0007  R55 
0008  R55 
.  .  0009  -.RSS 

t  SRZ  ,  0010  R55 

:  001 1  R55 

. .  0012  R55 

0013  R55 
0014  R55 
00151R55 
00 1 52R55 
001 53R55 

.  0016  R55 

0017  R55 
0018  R55 

. .  , ,  .  .  ,0019.  R55 

0020  R55 
0021  R55 

. . .  . . 002.2.  .R.5.5 

0023  R55 
0024  R55 

. .  .  ..Q025...R55 

0026  R55 
0027  R55 

. . .  00,28  R55 

0029  R55 
0030  R55 

. . . .  Q.031..R55 

0032  R55 
0033  R55 

_  . 0.0.34.  R5.5. 

0035  R55 
0036  R55 

. . . .  ..  ..003.7„  R55 

0038~  R5S 
0039  R55 

_ _  _ 0040  R5.5 

0041  R55 
0042  R55 

_ _  _  0043.  R.55 

0044  R55 

PRE  NOT  PUPILABLE90045  R55 

- _ D.Q16...B.5.5 

0047  R55 
0Cr48  R55 

_ 0049, R5.5. 

0050  R55 
0051  .R55 

.  . .  .  00.52  .R55 

0053  .R.55 
0054  R55 

_  _  __0.0.5,5.^5,5 
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FORTRAN  STATEMENTS  (cont'd.) 


900  IF  C  N0P1  +  N0P2  5  1000.  1100.  1200 

1000  PRINT  1001 

1001  FORMfiT  C37H  ERROR  ON  SECOND  CfiRD  OF  PROBLEM  DECK? 

CeLL  END9  . . . 

1  100  NM  =  320 

GO  TO  2000 

1200  NM  =576'  '  ....  .... 

CONSTRUCT  fl  LARGER  MATRIX  FROM  SCI.J? 

2000  DO  2022  1=1. NM  . 

DO  2022  J=1 ,8 

IF  CI-8?  2010,  2010,  2012 

2010  SSCI.J?  =  SCI.J? 

GO  TO  2022 

2012  SSCI.J?  =  SS  CI-8,  J? 

2022  CONTINUE 
65  FORMAT  ':8E10.4? 

IF  C  SENSE  SWITCH  2  ?  2023,  2024 

2023  PRINT  65,  CC  SSCI.J?,  J=1,S  ?,I=1,NM  ?  ... 

202h  do  2027  I=1,NM 

DO  2027  J=1.8 

.SSCI.J?  =  6.2831853  ♦  SS.CI.J?  . . 

2027  CONTINUE 

IF  CSENSE  SWITCH  2  ?  2028.  2029 

2028  PRINT  65,  CCSS Cl ,  J  ?  , J  =  1, ,.8  ?,  1  =  1,  NM  i . 

2029  IF  CN0P1  )  1000.  2150,  2030 

COMPUTE  ANGULAR  DISTRIBUTION.  OF  FLUK,..„.  _ 

2030  DO  2035  1=  1,576 

DO  2035  J=  1,8 

fad  CI,J?  =  SSCI.J?  *  FLUXA  Cl,.. I? 

FAD  ci,j:-  =  FA0CI,J?  *  ACJ? 

2035  CONTINUE 

C  SUMMATION  .  . . .  _ 

DO  2040  I=i,576 

FOCI?  =  0.0 

DO  2040  J  =  1 ,8  ...  _ _ _ 

. FOCI?  =  FDci?  +  fad  ci,j:) 

2040  continue 

,Dp,205.Q_.I  =  1  ,  18. _  _ _ 

DO  2050  j=1 , 4 

DO  2050  K=1,8 

_L  =_K _ +_  8.«_.  CJr  0_...+..  32,  *..CI.-.I? _ 

SOS  ri.j.KJ  =  FDCL? 

2050  CONTINUE 

C  PRINT  OUT  RESULTS  AS  4  18X8  MATRICES. _  _ 

DO  2100  J=i , 4 
PRINT  2060 

2060  FORMAT  C29H1  ANGULAR  DI_STRIBUT10N  OF  FJ,UX? _ 

'  PRINT  2070  “,  J 

2070  FORMAT  C22H0SEPARAT I  ON  DISTANCE  A, 11? 

PRINT  2080  _  .  _ 

2080  format  CIOOHO  K  EoT  E02  E03 

1  EOS  E06  E07  EOS  ? 

_D.Q  2100  1  =  1,18  _  _  _  _ _ 

PRINT  2*090  ,  Cl.  CSOS  CI.J.K?,  K=1  , 8  ?? 

2090  FORMAT  C1H  ,I2,8E12.40 

2.100  J:6nTINUE  _  _  _ _ _ 


. 00.56  R55, 

0057  R55 
0058  R55 

.  . . 0.059  .R  5.5 

0060  R55 
.  0061  R55 

_  .....  .00.62 -R55 

0063  R55 
0064  R55 

_  _ 0.065.  i{5.5 

0066  R55 
0067  R55 

. . .  ..  ..  .0.06.8.  R55 

0069  R55 
00701 R55 

. 0.0702.R55 

00703R55 

0071,1R55 

. aO,7.t.2R5.5. 

0.0.7J3R55 
.  .  .Q0714E55 

_ a0.1ZLB55 

.  ..0il.722R55 
.  ...D.Q723R55 

_ _ CLQ.i24fi.5.5 

.  ,  . 00.73  .R55 

Q074  R55 

_ '  0075  R55 

0076  R55 
0077  R55 

_ £LQ.7.8_E55 

,0079  R55 
0080  R55 

_ _ 0-08.1..  R5.5 

0082  R55 
0083  R55 

_ 0084  R5s 

0085  R55 
0086  R5b 

_ ! _ 0.0.87_.R5.5 

0088  R55 
0089  R55 

_ ^Q.CL9Q.„B.55 

0091  'r5S 
0092  R55 

_  00.93  R55 

0094  R55 
0095  R55 

_ Q096  R5;j 

0097  ’R55 
0098  R55 

_ -..JiCi3.2_R.55 

E04.  OiOO  R55 

0101  R55 

- - 0..lQ2.Ji.55. 

0103  R55 
.  01 64  R55 

_ 1 _ ^QlOSjai 


FORTRAN  STATEMENTS  (cont'd.) 


2150  IF  CNOP  2  >  1000,  3000,  2200 

C 

C  COMPUTE  flNQULflP  DISTRIBUTION  OF  DOSE  RATE. 

2200  DO  22.1.0  . .  ,  I ...=  .  1, 576 
DO  2210  J  =  1  ,8 

FflO  CI,J5  =  SSCI.J5  *  DCSEfl  CI,J> 

FflD  Cl.  J?  .  =  .  FHDCI,  J:)  * 

2210  continue 
c  slimmrtion 

DO  2220  1=1 ,576  . 

FDCn  =  0.0 

DO  2220  J=1,8 

FOCI?  =  FDCn  +  FRD  CI.jD 

2220  continue 

DO  2230  1=1,18 

DO  2230  J=1 ,4 

DO  2230  K:=1,8 

L  =  K  +  8  *  CJ-i)  +  32  *  CI-1> 

SOSCI.J.K:)  =  PDCLO 
2230  CONTINUE 

C  PRINT  OUT  RNGULRR  DISTRIBUTION  OF  DOSE  RPTE 

DO  2250  J=1,4 

PRINT  2240 

2240  FORMAT  C38H 1 RNGULfiR  DISTRIBUTION  OF  THE  DOSE  RRTE> 
PRINT  2.070,  J  . 

PRINT  2080 
DO  2250  1=1,18 

PRINT  2080  ,  Cl  ,  CSOS  CI,J,K;>,  K=1,S;:'  ? 

2250  CONTINUE 

IF  CNOP  3  >  1000,  50,  3000 

C 

C  COMPUTE  ENERGY  SPECTRUM  OF  FLUX 

3000  DO  3035  1=1,320 

DO  3035  J=1,8 

SSCI,JO  =  SSCI,J?  *  FLUXE  CI,.J? 

SSCI,J5  =  SSCI.J?  *  RCJ? 

3035  CONTINUE  ,  ,  ,,  . . . 

C  SUMMfiTION 

DO  3040  1=1 ,320 

FDCn  =  0.0  . . . 

DO  3040  J=1 ,8 

FDCn  =  FDCO  +  SS  CI,J8 

3040  continue  . .  .  ...  .  . 

DO  3050  1=  1,10 

DO  3050  J.=  1,4 

DO  3050,  K=  1,,8  •  „  . . . . 

L  =  K  +  8  ♦  CJ-»1J  32  *  C-I-O 

SOS  CI,J,K  >  s  FD<L> 

3050  continue  _ _ _ _ _ 

C  PRINT  OUT  ENERGY  SPECTRUM  OF  FLUX 

DO  3070  J=1 ,4 

PRINT  3060  .  _ _ 

3060  FORMAT  C28H1 TOTAL  FLUX  IN  ENERGY  GROUP  K> 

C 

PRINT  2070  ,  J  „  _  . 

PRINT '2686 
DO  3070  1=1,10 

PRINT  2090  ,  C  J, _ CSC5CI,U,K:),K=r,8i.  5  _ 
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0116  R55 
0117  R55 
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0119  R55 
0120  R55 
0121  R55 
0122  R55 
0123  R55 
...  .0124  R55 

0125  R55 
0126  R55 
0127  R55 
0128  R55 
0129  R55 
.  0130  .R55 
0131  R55 
0132  R55 
.  0133. R55 
0134  R55 
0135  R55 
01.36  R.55 
0137  R55 
0138  R55 
.  .  01,3,9  R  5.5 
0140  R55 
0141  R55 
.  ,0142  R55 

0143  R55 
0144  R55 
. .  0145  R55 
0146  R55 
0147  R5b 
.  .0148  R55 
0149  R55 
0150  R55 
.  0151  R55 
0152  R55 
0153  R55 
_0154  R55 
0155  R55 
0156  R55 
0157  R55 
0158  R55 
0159  R55 
g160  R55 
0161  R55 
.0162  R55 
.  __  0J63  R55 
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FORTRAN  STATEMENTS  (cont'd.) 


3070  CONTINUE 
4000  GO  TO  50 
ENDC2. O 


0164  R55 
0165  k55 
0166  R55 
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